The three-flavor chiral phase structure in hot and dense QCD matter 
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Chiral symmetry restoration at nonzero temperature and quark densities are investigated in the 
framework of a linear sigma model with Nf = 3 light quark flavors. After the derivation of the grand 
potential in mean-field approximation, the nonstrange and strange condensates, the in-medium 
masses of the scalar and pseudoscalar nonets are analyzed in hot and dense medium. The influence 
of the axial anomaly on the nonet masses and the isoscalar mixings on the pseudoscalar rj-ri' and 
scalar (t(600)-/o(1370) complex are examined. The sensitivity of the chiral phase transition as well 
as the existence and location of a critical end point in the phase diagram on the value of the sigma 
mass is explored. The chiral critical surface with and without the influence of the axial U(1)a 
anomaly is elaborated as a function of the pion and kaon masses for several values of the sigma 
mass. 

PACS numbers: 12.38.Aw, 11.30.Rd,14.40.Aq 



I. INTRODUCTION 

The understanding of the properties of strongly inter- 
acting matter under extreme conditions is one of the most 
fascinating and challenging tasks. General features of hot 
and dense matter are summarized in the QCD phase dia- 
gram which can be probed by ultrarelativistic heavy ion 
experiments such as the RHIC (BNL), LHC (CERN) and 
the planned future CBM experiment at the FAIR facility 
in Darmstadt. 

Theoretical considerations indicate that at high tem- 
perature and high baryon densities there should be a 
phase transition from ordinary hadronic matter to a chi- 
rally symmetric plasma of quarks and gluons Sev- 
eral issues concerning this transition are not yet clarified 
0. QCD in this temperature and density regime is a 
strongly coupled theory and hence perturbation theory 
cannot be used. In the absence of a systematically im- 
provable and converging method to approach QCD at 
finite density one often turns to model investigations sec 
e.g. 0]. These models incorporate the important chiral 
symmetry breaking mechanism of QCD but neglect any 
effects of confinement. Only recently, certain aspects of 
confinement based on the Polyakov loop have been incor- 
orated in chiral effective models in a systematic fashion 
|5|| and interesting conclusions could be drawn (sec 
e.g. i). 

The most prominent finding from low-energy chiral ef- 
fective models is the QCD critical end point (CEP) Q. 
Common to almost all effective model calculations is that 
the chiral phase transition is continuous in the low den- 
sity region and discontinuous in the high density regime. 
Consequently, the endpoint where the phase transition 
ceases to be discontinuous is the QCD critical end point. 
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Unfortunately, several obvious and related features such 
as the exact location of this point in the QCD phase di- 
agram cannot be predicted by these models. 

On the other side, lattice QCD simulations are im- 
portant alternatives to effective models calculations and 
can gain much insi ghts in the QCD phase structure 
B i,M HH, E, El, in, Due to the notorious sign 
problem emerging at finite baryon density reliable predic- 
tions for QCD are still not conclusive. Even worse, re- 
cently different lattice methods that circumvent the sign 
problem are in confiict to each other. For example, us- 
ing the imaginary chemical potential method for three 
physical quark masses no critical endpoint in the phase 
diagram is found [H, [13, M, M,M,^- 

The present work is an extension of a previous analysis 
within an effective linear sigma model (LcrM) with two 
quark fiavors to three quark ffavors [2^ . [23j . The restora- 
tion of the chiral SU{3) x SU{3) and axial U{1)a sym- 
metries with temperature and quark chemical potentials 
are investigated. The axial U{1)a anomaly is considered 
via an effective 't Hooft determinant in the Lagrangian 
which breaks the U{1)a symmetry. The restoration of 
the U {1)a symmetry is linked to a vanishing of the topo- 
logical susceptibility which can further be related to the 
77' mass via the Witten-Veneziano relation [U, [2^ . 

Some results depend sensitively on the model param- 
eters which are tuned to reproduce the vacuum phe- 
nomenology. There are model-input parameters such as 
the (T meson mass which are poorly known experimen- 
tally. The generic findings of several parameter fits over a 
broad range of input parameters are compared. Further- 
more, the extrapolation towards the chiral limit is also 
addressed and the mass sensitivity of the chiral phase 
structure is investigated. 

The paper is organized as follows: after introducing the 
LcrM with three quark flavors, some symmetry breaking 
patterns in the vacuum are briefly discussed. In Sec. IIIII 
the grand thermodynamic potential is derived in mean- 
field approximation. In Scc. lIVI a discussion of the model 
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parameter fits is given. Since the experimental situation 
concerning the scalar cr-meson, ct(600), is not settled we 
consider a wide range of different values of the cr-mcson 
mass, m^, as input parameter. All parameter sets are 
collected in App. [X) 

For the LcrM without quarks it is known that the stan- 
dard loop expansion and related approximation methods 
at finite temperatures fail and imaginary meson masses 
are generated. In our approximation no such artifacts 
occur which enables us to perform a careful and detailed 
analysis of chiral symmetry restoration in hot and dense 
matter. This is demonstrated in Sec. |V] where the pseu- 
doscalar and scalar meson masses at finite temperatures 
and chemical potentials with and without axial U{1)a 
symmetry breaking are investigated. All mass expres- 
sions are summarized in App. |B1 In addition, the scalar 
and pseudoscalar flavor mixing behavior in the medium 
is explored. Various definitions are deferred to App. [C] 

The grand potential determines all thermodynamic 
properties. The resulting phase diagrams are presented 
in Section IVII where the mass sensitivity of the chiral 
phase boundaries is also explored. Subsequently, the 
shape of the chiral critical surface which confines the re- 
gion of the chiral first-order transitions in the mT^-mK 
plane at the critical chemical potential is evaluated for 
several values of ttIo-. Finally, in Sec. IVIII a summary 
with concluding remarks is given. 

II. LINEAR SIGMA MODEL WITH THREE 
QUARK FLAVORS 

The Lagrangian, Lqm ~ Cq + C„i, of the SU{3)l x 
SU{3) R symmetric LctM with three quark flavors consists 
of the fcrmionic part 

Cq = q {i^ - g Ta {aa + ijsT^a)) q (1) 

with a flavor-blind Yukawa coupling g of the quarks to 
the mesons and the purely mesonic contribution 

Cm = Tr - m2 Tri^^) - Ai [Tr(0t0)] ' 

Tr (0^0)^ + c (det((/)) + det((/)^)) 
+ Tr + 0t)] . (2) 

The column vector q = (u, d, s) denotes the quark field 
for Nf = 3 flavors and A^c = 3 color degrees of freedom 
[2^. The 0- field represents a complex (3 x 3)-matrix and 
is defined in terms of the scalar (Tq and the pseudoscalar 
tTq meson nonet 

(f> = Ta(j)a = Ta {(Ja + ilTa) ■ (3) 

The Ta = Aa/2 with a = 0, . . . , 8 are the nine generators 
of the /7(3) symmetry, where the Aa are the usual eight 

Gell-Mann matrices and Aq = 1. The generators 
Ta are normalized to Tr(TaTb) ~ Sab/'^ and obey the 
C/(3) algebra [Ta,Tb] = ifabcTc and {Ta,Tb} = dabcTc 



respectively with the corresponding standard symmetric 
dabc and antisymmetric fabc structure constants of the 
SU{3) group and 




Chiral symmetry is broken explicitly by the last term in 
Eq. ^ where 

H = Taha (5) 

is a (3 X 3)-matrix with nine external parameters ha- In 
general, one could add further explicit symmetry break- 
ing terms to Lm which are non-linear in (p [2?! . [28j but 
this is ignored in this work. 

Due to spontaneous chiral symmetry breaking in the 
vacuum a finite vacuum expectation value of the </> field, 
(j), is generated which must carry the quantum numbers 
of the vacuum [2§|. As a consequence, only the diago- 
nal components hojh^ and /ig of the explicit symmetry 
breaking term can be nonzero. This in turn involves three 
finite condensates (Tq, CT3 and as of which (T3 breaks the 
SU{2) isospin symmetry. In the following we shall re- 
strict ourselves to a 2 -|- 1 fiavor symmetry breaking pat- 
tern and neglect the violation of the isospin symmetry. 
This is refiected by the choice ho 7^ 0, /13 = 0, /is 7^ and 
corresponds to two degenerated light quark flavors (u, d) 
and one heavier quark flavor (s). 

Besides the explicit symmetry breaking terms /iq and 
/ig the model has five more parameters: the squared tree- 
level mass of the meson fields m^, two possible quar- 
tic coupling constants Ai and A2, a Yukawa coupling g 
and a cubic coupling constant c which models the axial 
U{1)a anomaly of the QCD vacuum. The U{1)a symme- 
try of the classical QCD Lagrangian is anomalous [sot . 
i.e. broken by quantum effects. Without the anomaly 
a ninth pseudoscalar Goldstone boson corresponding to 
the spontaneous breaking of the chiral U{3)l x U{3)ii 
symmetry should emerge. However, experimentally, the 
lightest candidate for this boson is the 77' meson, whose 
mass is of the order m,f ^ 960 MeV which is far from 
being a light Goldstone boson. The explicit breaking of 
the U{1)a symmetry is held responsible for the fact that 
the Tj' mass is considerably larger than all other pseu- 
doscalar meson masses. This well-known U{1)a problem 
of QCD is effectively controlled by the anomaly term c 
in the Lagrangian. The comprehensive procedure of how 
to fix the parameters will be given in Sec. IIVI 

Depending on the signs and values of the parame- 
ters several possible symmetry breaking patterns in the 
vacuum can be obtained (see also [28[ for more de- 
tails). Without explicit symmetry breaking terms, i.e. for 
H ^ 0, and without an explicit U{1)a symmetry break- 
ing term, i.e. for c — 0, the Lagrangian has a global 
SU{3)v X Ui3)A ^ SU{3)v x SU{3)a x U{1)a symme- 
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try if the quartic coupling A2 and are positive^. If 
the mass parameter m? changes sign the symmetry is 
spontaneously broken down to SU{3)v- The first quar- 
tic coupling Ai has no influence on the symmetry break- 
ing. Because of the breaking of the U{3)a symmetry 
nine pseudoscalar Goldstone bosons arise which form the 
entire nonet consisting of three pions, four kaons, the 77 
and ?/ meson. The scalar nonet belongs to the SU{3)v 
group which has a singlet and an octet representation. 
All masses of the octet particles are degenerate. The a 
meson belongs to the singlet and its mass is in general 
different from the masses of the other octet particles. 

By setting c 7^ the effects of the U{1)a symmetry 
breaking, caused by a nonvanishing topological suscepti- 
bility, are included and the symmetry of the Lagrangian 
is reduced to SU{3)v x SU{3)a- Due to the spontaneous 
symmetry breaking of the SU{3)a, the vacuum has a 
SU{3)v symmetry. In this case the entire pseudoscalar 
octet is degenerated and only eight Goldstone bosons ap- 
pear. The 77' meson, the would-be Goldstone boson, is 
still massive in this case. The masses of the scalar parti- 
cles are not modified by the U{1)a breaking. 

With explicit symmetry breaking terms, more pre- 
cisely, for only finite ho and hg terms, the vacuum 
SU{3)v symmetry is explicitly broken down to the 
isospin SU{2)v symmetry since the /13 term is set to zero. 
This symmetry pattern is already a good approximation 
to nature because the violation of the isospin symmetry 
is small anyway. The resulting ground state spectrum for 
this symmetry pattern will be discussed in Sec. IV Bl 



III. GRAND POTENTIAL 

In this section, the derivation of the grand thermody- 
namic potential for the three-flavor model is given. We 
will use a mean-field approximation similar to the one 
for the two-flavor model in . The mean-field approxi- 
mation is simple in its application, in particular at finite 
temperature and quark densities. Low-energy theorems, 
such as e.g. the Goldstone theorem or the Ward identities 
arc also fulfilled at finite temperatures and densities. In 
this way we can circumvent more advanced many-body 
resummation techniques which are usually necessary to 
cure the breakdown of naive perturbation theory due to 
infrared divergences. For example, it is well-known that 
the standard loop expansion or related expansion meth- 
ods of the SU{3) version of the LcrM with or without 
quarks break down at finite temperature and imaginary 
meson masses are generated in the spontaneously bro- 
ken phase [m, m, |33,, i35(] . Contributions of thermal 



More precisely, without the determinant the Lagrang ian i^qm is 
(7(3) X f/ (3) invariant which is isomorphic to S';7(3)i^ X 5f7(3) n X 
U{1)b X (7(1) A- The U{1)b is related to the baryon number and 
always conserved which is why we have neglected it. 



excitations to the meson masses are neglected in these ap- 
proximation schemes which result in a too rapid decrease 
of the meson masses and the e.g. squared pion mass be- 
comes negative for temperatures much below the phase 
transition. This deficiency can be cured by self-consistent 
resummation schemes such as e. g. t he Hartree approxi- 
mation in the C JT formalism [s^jUBl or the so-called Op- 
timized Perturbation Theory (OPT) e.g. [13] and variants 
thereof. Recently, the OPT method has been frequently 
applied to the three flavor LcrM with and without quarks 
at finite temperature and baryon densities [ssl . [39l | . How- 
ever, the predictive power of the OPT method depends 
on how it is implemented and approximations thereof 
are made. For instance, when the external momentum of 
the self-energy are taken on-shell a solution of the corre- 
sponding gap equation and also of the equation of state 
cease to exist above a certain temperature, particularly 
below the critical one. Details and some improvements 
of certain appro ximations in the OPT framework can be 
found in [SSl l39l|. 

All these problems do not emerge in the mean-field 
approximation used here. This enables us to study the 
phase structure of the more involved three flavor model 
in great detail and in a rather simple framework. 

In order to calculate the grand potential in mean-fleld 
approximation we start from the partition function. In 
thermal equilibrium, the grand partition function is de- 
fined by a path integral over the quark/ antiquark and 
meson fields 



Z^JY[ '^'^a'DTTaJ VqVq exp f - y drj d^xC 



(6) 



where T is the temperature and V the three-dimensional 
volume of the system^. For three quark flavors the Eu- 
clidean Lagrangian generally contains three indepen- 
dent quark chemical potentials /i j 



Due to the assumed SU{2)y isospin symmetry we neglect 
the slight mass difference between an u- and d-quark and 
the light quark chemical potentials become equal. In the 
following we denote the degenerated light quark quanti- 
ties by an index q, i.e. the light quark chemical potential 
by = flu = fJ.d, and the strange quark quantities by an 
index s. 

The calculation of the partition function in the mean- 
field approximation is performed similar to Refs. p^ . 
[ioj for the two-flavor case. The quantum and ther- 
mal fluctuations of the mesons are neglected and the 
quarks/ antiquarks are retained as quantum fields. This 
means that the integration over the mesonic fields in 



An irrelevant normalization constant is suppressed. 
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Eq. ([6]) is dropped and the fields are replaced by their 
non-vanishing vacuum expectation values (f> = Tq^o + 
TgCTs. The remaining integration over the Grassmann 
fields yields a determinant which can be rewritten as a 
trace over a logarithm. Evaluating the trace within the 
Matsubara formalism, the quark contribution Qqq{T, fij) 
of the grand potential is obtained [4l| . The ultraviolet 
divergent vacuum contribution to flqq{T, fif) which re- 
sults from the negative ei ierg y states of the Dirac sea has 
been neglected here (cf. [22|, [40| for further details) . Fi- 
nally, the total grand potential is obtained as a sum of 
the quark contribution and meson contribution [/((Tq, (Tg) 
as 



-TlnZ 
V 



Uiao,as) + nqqiT,tif) . (7) 



Explicitly, the quark contribution reads 



oo 



f — U,d,S ) 



(Pk 
f27r)3 



+ Hl-nq,J{T,^if))} (8) 



with the usual fermionic occupation numbers for the 
quarks 



The meson potential modifies accordingly 



U{(Tx,(7y) 



{"'I + - h.j,ax - hyffy 



+ ^<y'H + \ (2Ai + A2) <jt + \ (2Ai + 2A2) <jI , (13) 



therein the explicit symmetry breaking parameters 
and /ig have also been transformed according to Eq. 
The order parameters for the chiral phase transition are 
identified here with the expectation value a^. for the non- 
strange and with dy for the strange sector. They are ob- 
tained by minimizing the total thermodynamic potential 
in the non-strange and strange directions 



on 

dax 



day 



(14) 



The solutions of these coupled equations determine the 
behavior of the chiral order parameters as a function of 
T and chemical potentials, fiq and fis- Note, that the 
in-medium condensates are also labeled with a bar over 
the corresponding fields. 



+ exp((i?,j-M/)/r) 



(9) 



and antiquarks nqj{T, fif) = nqj{T, —fif) respectively. 
The number of internal quark degrees of freedom is de- 
noted by i^c = 2iVt, = 6. The flavor-dependent single- 
particle energies are 



EqJ = Jk 



(10) 



with the fiavor-dcpendcnt quark masses mf which are 
also functions of the expectation values ctq and CTg . 

The vacuum condensates (Tq and ag are members of 
the scalar [J^ = 0+) nonet and both contain strange 
and non-strange components. For the further analysis 
it is more convenient to convert the condensates into a 
pure non-strange and strange part. This is achieved by an 
orthogonal basis transformation from the original octet- 
singlet basis (ctq, cTg) to the non-strange (era;) and strange 
(CTy) quark fiavor basis 



V3\i -A 



CTO 



(11) 



As a consequence, the light quark sector decouples 
from the strange quark sector (cf. e.g. [3^) and the quark 
masses simplify in this new basis to 



Too 



: gax/2 



rus 



gay 



/V2 



(12) 



IV. PARAMETER FITS 



The LaM with three quark fiavors has altogether seven 
parameters m^, Ai, A2, c, g, Hq, /ig and two unknown 
condensates ax and ay. The six parameters m^, Ai, A2, 
c, hx and hy of the mesonic potential are fixed in the 
vacuum by six experimentally known quantities. Simi- 
lar to Ref. [13], we have chosen as input the low-lying 
pseudoscalar mass spectrum, m^r and rriK, the average 
squared mass of the 77 and 77' mesons, tti^ -I- rn^, , and 
the decay constants of the pion and kaon, /^r and /k, 
and in addition the scalar a meson mass 771^. We can 
then predict the scalar meson masses mag, tti^, mj^, the 
difference of the 77,77' squared masses, rn^^ — 
scalar and pseudoscalar mixing angles 6s, Sp- 
in analogy to e.g. Ref. [s^l the values of the conden- 
sates are determined from the pion and kaon decay con- 
stants by means of the partially conserved axial-vector 
current relation (PC AC). In the strange-non-strange ba- 
sis they are given by 



777'^/ and the 



(15) 



The average squared 77 and 77' meson mass determines 
the parameter A2 by 
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3(2 fK ' UWk - C^fK + U)ml - 1{m%, + m^)(/K - h) 

{3f^ + 8fKifK-f.))ifK-frr) 



(16) 



r 



and the U{1)a anomaly breaking term c is fixed by A2 
and the diff'erence of the pion and kaon masses squared 
via 



— — M{2fK — Jtt) 



(17) 



Note, that without anomaly breaking, i.e. for c = 0, the 
average ry-ry' meson mass is not used anymore for fixing 
the parameter A2. It is then given by the kaon and pion 
masses and decay constants only, 



(18) 



The input parameters from the pseudoscalar sector in- 
volve only a relation between Ai and vn? . Therefore, fur- 
ther input from the scalar sector is necessary. In principle 
two possible options are available. At first, the parame- 
ter is expressed as a function of the yet undetermined 
parameter Ai. This can be achieved by using the equa- 
tion for the pion mass or the kaon mass (sec App. IB|) . 
In this way the dependence of the a meson mass 
(or of the /o(1370) meson mass) can be transformed in 
a Ai dependence since the scalar mixing angle does not 
depend on w? . By fixing the mass of the a meson (or 
of the /o(1370) meson) Ai is determined by solving the 
corresponding equation. Afterwards, the parameter 
follows immediately, since Ai is fixed. 

The explicit symmetry breaking terms and hy in 
the non-strange-strange basis are related to the pion and 
kaon masses by the Ward identities 



hy = V2fKm\ 



V2 



(19) 



These relations can be derived by using the gap equa- 
tions, Eqs. (dH). The last open parameter, the value of 
the Yukawa coupling is fixed from the non-strange con- 
stituent quark mass 



g = 2mq/cFx 



(20) 



For example, using for the light constituent quark mass a 
value of ruq = 300 MeV we obtain g ^ 6.5 and can predict 
a strange constituent quark mass nis k, 433 MeV. 

Since the experimental situation concerning the broad 
a (or /o(600)) resonance is not yet clear, cf. [i^, we 
will use different input values for rria in the range of 
rricr — 400 — 1000 MeV and will investigate its mass de- 
pendence on various quantities (see also [11]). In App. \K\ 
several parameter sets for different TOo- values with and 



without effects of the axial U{V)a anomaly are summa- 
rized fTab. [TT)l. Furthermore, a discussion of the param- 
eter sets with respect to spontaneous symmetry breaking 
can be found in this appendix. The corresponding pre- 
dictions of the scalar and pseudoscalar meson masses and 
mixing angles are collected in Tab. U 



V. CHIRAL SYMMETRY RESTORATION 

Having fixed the model parameters we can now evalu- 
ate the grand potential numerically. In the following we 
present our results for the chiral symmetry restoration 
at finite temperature and finite quark density with and 
without axial anomaly breaking. Throughout this sec- 
tion, the axial anomaly breaking term is kept constant, 
in particular, independent of the temperature and the 
chemical potentials. 



A. Condensates 

The chiral phase structure of the underlying three- 
fiavor model is completely governed by the total ther- 
modynamic potential. Hence, the solution of the gap 
equations determines the behavior of the conden- 
sates as a function of temperature and quark chemical 
potentials. In general, the three quark chemical poten- 
tials are independent but here we will consider symmet- 
ric quark matter and define a uniform chemical potential 

In Fig. [T] the nonstrange ax and strange ay conden- 
sates are shown as a function of temperature for van- 
ishing chemical potential /i for ra^r — 800 MeV. The 
reason for choosing this value for m„ will become clear 
later on, see also App. [X] The solid lines in this fig- 
ure are obtained with an explicit axial U{1)a symmetry 
breaking term while the dashed line corresponds to the 
anomaly free case, i.e. c = 0. The difference in the 
nonstrange condensate dx , caused by the anomaly, is not 
visible in the figure. The condensates start at T = 
with the fitted values, ax = 92.4 MeV and ay ~ 94.5 
MeV. The temperature behavior of both condensates 
shows a smooth crossover. The temperature derivative of 
the nonstrange condensate peaks around T ^ 181 MeV. 
The precise value of this pseudo-critical temperature de- 
pends on the value of m„ in the vacuum. For smaller 
values of m„ the pseudocritical temperature decreases 
(cf. Sec. IVI|) . The chiral transition in the strange sector 
is much smoother due to the larger constituent strange 
quark mass, nis = 433 MeV. As a consequence, the chiral 
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FIG. 1: The non-strange, ax and strange, ay condensates as a 
function of temperature for vanishing chemical potentials with 
(solid) and without U{1)a anomaly (dashed). The anomaly 
does not modify the non-strange condensate. 



can be seen in Fig. [2l At fic the phase transition is of 
first-order and three solutions of each gap equation pi]) 
appear corresponding to two degenerate minima and one 
maximum of the effective potential. Due to the explicit 
symmetry breaking both condensates remain always fi- 
nite in the symmetric phase. The phase transition is 
mainly driven by the nonstrangc condensate while the 
jump in the strange condensate is negligible. Above the 
transition for fi > fic and below the strange Fermi sur- 
face at /i = nis ~ 433 MeV the strange condensate stays 
constant. The axial anomaly has almost no influence up 
to the strange quark Fermi surface. Only for chemical 
potentials larger than /i ^ 433 MeV the strange con- 
densate melts faster if the U{1)a symmetry breaking is 
taken into account. For large chemical potentials this dif- 
ference vanishes again and both strange condensates will 
become identical. Furthermore, the in-medium behavior 
of the nonstrangc condensate is not modified by the 
anomaly. 
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FIG. 2; Similar to Fig[l]but as a function of /i for T — 0. 



SU{2) X SU{2) symmetry is restored more rapidly. With 
axial anomaly the strange condensate melts a little earlier 
but only for temperatures above the transition. At very 
high temperatures both condensates become degenerate, 
indicating chiral SU{3) x SU{3) symmetry restoration. 

If one uses a temperature dependent anomaly term by 
making use of lattice results for the topological suscep- 
tibility which yields e.g. a decreasing anomaly term for 
increasing temperatures, a faster effective restoration of 
the axial symmetry can be achieved, see e.g. [i^. |45|. 

For zero temperature and finite chemical potential 
both condensates are independent of fi in the broken 
phase until the Fermi surface of the light quarks is 
reached. For zero temperature the light Fermi surface 
coincides with the light quark mass, i.e. jj, = niq = 300 
MeV. Before the chiral transition takes place at a critical 
chemical potential ~ 352 MeV the nonstrangc con- 
densate drops by about 10% from its vacuum value as 



B. The scalar-pseudoscalar meson spectrum 

In the following the in-medium scalar and pseudoscalar 
meson mass spectrum is analyzed. The derivation of the 
in-medium masses as well as the mass formulae are col- 
lected in App. |B] 

We start with the discussion of the mass spectrum at 
nonzero temperature and vanishing quark chemical po- 
tential. The meson masses as a function of the tempera- 
ture for = are shown in Figs[3]and[H In the left panels 
of the respective figures the U{1)a symmetry breaking is 
explicitly taken into account while in the right panels the 
breaking is neglected. 

The masses of the pion and the a meson and also the 
masses of the if and the oq meson degenerate approx- 
imately at the same temperature Tc ~ 181 McV. This 
temperature behavior signals the effective restoration of 
chiral SU{2) x SU{2) symmetry in the non-strange sec- 
tor via a smooth crossover transition. The chiral partners 
(K, k) show a similar temperature behavior but degener- 
ate with the rj meson at a higher temperature T ^ 240 
MeV. At the chiral transition T ~ 181 MeV the k meson 
becomes lighter than the oq meson. In contrast to [131 
the rj' is always heavier than the kaon for all tempera- 
tures. Only the /o(1370) meson mass does not show a 
tendency to converge to the rj mass in the temperature 
region shown because chiral symmetry in the strange sec- 
tor is restored only very slowly. The intersection point of 
the fo and the r/ mass coincides with the inflection point 
of the strange condensate. Nevertheless, the fo meson 
will degenerate with the remaining meson octet at very 
large temperatures. 

The mass gap in the restored phase for T > Tc between 
the two sets of the chiral partners, (ct, tt) and (oq, r]'), i.e. 
TOtt = rria- < uiao — rrijy is a consequence of the U{1)a 
breaking term. This gap is generated by an opposite sign 
of the anomaly term in the scalar and pseudoscalar me- 
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FIG. 3: In-medium meson masses, (tt, a) and (77', ao), as a function of temperature for /i = with C/(1)a anomaly breaking 
(left panel). Without anomaly breaking (right panel) the r\ meson degenerates with the pion mass. See text for further details. 
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FIG. 4; Similar to Fig.|3]for the chiral partners (r\, fo) and {K, k) (left panel with U{1)a anomaly breaking). Without anomaly 
breaking (right panel) the 77 meson increases about 100 MeV in the vacuum. 



son masses, cf. App. [BJ It is basically given by \/2cay, 
i.e. proportional to the anomaly term c and the strange 
order parameter ay. The non-strange condensate ax is 
already negligible for temperatures above Tc- For higher 
temperatures (T ^ 400 MeV) the U{1)a symmetry gets 
effectively restored and the mass gap between the chiral 
partners will vanish. Finally, for very large temperatures 
compared to the strange quark mass the difference be- 
tween the strange and non-strange mesons disappear and 
all meson masses will degenerate. 

Without U{1)a symmetry breaking the mass gap be- 
tween the chiral partners, (cr, tt) and (ao,?/) vanishes in 
the restored phase and all four meson masses degenerate 
at the same critical temperature Tc ^ 181 MeV coin- 
ciding with the inflection point of the non-strange con- 
densate. Above this temperature the axial symmetry 
is restored but the full restoration of the C/(3) x J7(3) 
symmetry is still not yet completed because the chiral 
partners {K, k) degenerate with the 77 at a higher tem- 



perature T ^ 240 MeV. This temperature value and the 
value of Tc are not changed by the anomaly as expected 
since the non-strange condensate is not influenced by the 
anomaly. Interestingly, a recent mean-field study within 
the three-flavor NJL model with various effective U{1)a 
anomaly implementations found an explicit difference for 
the chiral non-strange transition temperatures with and 
without explicit U{1)a symmetry breaking (cf. Tab. Ill 
in El). 

As for the case with anomaly, the chiral partners (77, fo) 
degenerate but only for temperatures around 300 MeV 
because these mesons are purely strange states and chi- 
ral symmetry in the strange sector is very slowly restored. 
A mild anomaly dependence of the intersection point of 
the fo and the rj meson is observed. There is no inverse 
mass ordering of the 77 meson and the kaon at finite tem- 
perature as found in 3J| . In the vacuum the mass of the 
fo increases by about 60 MeV if the anomaly is neglected. 

Without the anomaly term the 77' meson degenerates 
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with the pion already in the vacuum and stays degener- 
ated with the pion for all temperatures. Hence, in the 
vacuum the mass of the t]' drops down considerably from 
963 MeV to 138 MeV. In fact, it has been shown that 
the mass of the v( must be less then \fivn.,^ ^ 240 MeV 
if the C/(1)a symmetry is not explicitly broken [30| . 

In general, one can summarize the mass spectrum in- 
medium behavior in the following way: the bosonic ther- 
mal contributions decrease the meson masses while the 
fermionic parts increase the masses. For small tem- 
peratures the quark contribution is negligible and for 
high temperatures it dominates the mesonic contribution 
yielding rising and degenerate meson masses. 

All meson masses are controlled by the two explicit 
symmetry breaking parameter /i^, and hy. They arc 
determined by the tree-level Ward identities, Eq. 
which guarantee the Goldstone theorem at zero temper- 
ature: for vanishing external parameter the pion mass 
must also vanishes because is then finite. In this case, 
the other symmetry breaking parameter hy generates 
only a finite value for the kaon mass. Furthermore, the 
chiral limit can be reached by setting all explicit symme- 
try breaking parameters to zero. But in order to obtain 
finite vacuum expectation values for the condensates the 
symmetry must be spontaneously broken. This requires a 
negative parameter. Later, we will use several param- 
eter fits for various values of the sigma meson mass which 
partly have a positive rr? parameter (see App. For 
these parameter sets one cannot reach the chiral limit by 
just setting the explicit symmetry breaking parameters 
to zero. But these parameter sets are still well suited for 
fitting the physical mass point. For instance, choosing a 
m„ — 400 MeV the parameter fit results with or without 
anomaly in a positive vr? parameter and the chiral limit 
cannot be reached for this parameter set. In [33 | another 
strategy to investigate the chiral limit was adopted by 
performing a separate extra parameter fit where an av- 
erage of the experimental mass values in the scalar octet 
spectrum together with some extrapolated quantities to- 
wards the chiral limit as input have been used. However, 
all in all the extrapolation towards the chiral limit re- 
mains questionable for both procedures. 

For the parameter set with e.g. rria = 800 MeV we can 
reach the chiral limit. For c ^ we obtain a massless 
pseudoscalar octet and a finite m^/ = 767 MeV due to the 
\J(X)a symmetry breaking. All scalar octet masses are 
degenerate at 840 MeV and the mass of the sigma meson 
drops down to 620 McV. Moreover, without \J(\)a sym- 
metry breaking all nine pseudoscalar mesons are massless 
and the scalar octet masses are degenerated into 780 MeV 
and m^r = 712 MeV. 

In our approximation the Goldstone's theorem is also 
valid at finite temperature and chemical potentials, 
meaning that in the chiral limit the masses of the Gold- 
stone bosons stay massless in the broken phase. Even in 
the presence of quarks both Ward identities in arc 
always fulfilled for all temperatures and quark chemical 
potentials which can be shown analytically. 



40 
30 
20 
10 

-10 
-20 
-30 
-40 
-50 
-60 



0pW/U(1)A 
0pW/oU(1)A 

03w/U(1)a 
©s w/o U(1)a 



50 100 150 200 250 300 350 400 
T [MeV] 



FIG. 5: The scalar Qs and pseudoscalar Qp mixing angles as 
a function of temperature for = without and with Ua(X) 
anomaly. 



Another observation in the Figs [3] and |4] is the tem- 
perature behavior of the scalar a and /o meson around 
T ^ 325 MeV. There is a kink visible in the curves and 
the meson masses seem to interchange their identities for 
higher temperatures. In order to elucidate this behavior 
we analyze the scalar and pseudoscalar mixing angles in 
the following. 



C. Flavor mixing at finite temperature 

The investigation of the mixing angles of the scalar and 
pseudoscalar isoscalar states provides further insights of 
the axial C/(1)a symmetry restoration. In order to clar- 
ify our findings some underlying definitions and relations 
between different bases are given in App. [Cl Both mix- 
ing angles, the pseudoscalar Qp and the scalar one 
are shown in Fig. [5] as a function of temperature for 
/i = with and without explicit U{X)a symmetry break- 
ing. In the broken phase, i.e. for T < 200 MeV a strong 
influence of the anomaly on the pseudoscalar sector is 
found while almost no effect is seen in the scalar sec- 
tor. With anomaly the nonstrange and strange quark 
states mix and generate an pseudoscalar mixing angle 
6'p ~ —5° at T = 0. For increasing temperatures the 
mixing angle stays almost constant in the chirally bro- 
ken phase. Around the chiral restoration temperature 
Tc ^ 180 McV the angle increases smoothly towards 
the ideal mixing angle Qp arctanl/V2 ~ -1-35° corre- 
sponding to (j)p = 90° , where (j)p denotes the pseudoscalar 
mixing angle in the strange-nonstrange basis (see App.lCl 
for details). At high temperatures this means that the 
?7 meson becomes a purely str ang e and the 77' a purely 
non-strange quark system (cf. [46[). 

This is also demonstrated in Fig. [S] where the physi- 
cal 77-77' and the nonstrangc-strange T^nS'Vs complex are 
shown as a function of temperature for = 0. At T = 
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the nonstrange mass of the 77NS meson, m,j„g , is larger 
than the strange mass 77i,,g since the pseudoscalar mixing 
angle is larger than 0^ = 45°, respectively Qp = —9.74°. 
For the mixing angle of Bp ^ —5° (0^ ^ 49.74°) we ob- 
tain m^pjs ~ 813 MeV and 777.^3 ~ 746 MeV. 

At the chiral transition temperature Tc ~ 180 MeV the 
77' meson becomes purely nonstrange (77' — > 77ns) while 
the 77 becomes a purely strange quark system (77 77s). 
In this temperature region the mixing angle grows to the 
ideal dp +35° (respectively 0p 90°). No crossing 
of the 77NS and 77s or anticrossing'^ of the physical 77-77' 
complex is observed for all temperatures since </>p(r) is 
always above 45° [i^. 

Without anomaly the pseudoscalar mixing angle is al- 
ready ideal for zero temperature and stays ideal for all 
temperatures, i.e. Bp ^ -1-35°. This means that already 
at r = the 77 and 77' mesons are ideally flavor-mixed. 
The 7;' is a purely light quark system and the 77 is a 
purely strange quark system. Without anomaly the 77' 
degenerates in mass with the pion. Hence, the 77' be- 
longs to the class of nonstrange particles. The order- 
ing of the corresponding nonstrange-strange masses is 
reversed, i.e. without anomaly is larger than vn^-^^^ 
since Tn^g = 771^ and ra^^^ = m^'. 

In the scalar sector the mixing angle 9s shows no in- 
fluence of the axial anomaly in the broken phase (see 
Fig. [5]). In both cases, with and without anomaly break- 
ing, the mixing angle is almost ideal 63 ^ +31° at T = 
but the precise vacuum value depends strongly on the 
value for the fitted scalar sigma meson mass in contrast 
to the pseudoscalar angle which is independent of m^, 
see Tab. U in App. [X] For increasing 777^ the scalar 



A crossing of the rjNs and r]s masses corresponds to an anticross- 
ing of the physical rj and r)' masses via Eqs. l|C9[l - l|C10[l (cf. also 



mixing angle 6*5 also increases at T = 0. As a con- 
sequence, for larger 77^0- the nonstrange-strange crNs-crs 
complex and the physical a and /o(1370) mesons degen- 
erate more and more at T = 0, meaning that the a me- 
son tends to a pure nonstrange quark system, a — !■ ctns, 
and the fo to a pure strange system, /o as- For in- 
stance, we obtain for ma{mfg) = 400(1257) MeV respec- 
tively TTicTNs ('^crs ) = 561(1131) MeV and for ttIo- (tti/q ) = 
800(1341) MeV m^^^im^^) = 804(1276) MeV. 

At the chiral transition Bs grows again to the ideal 
one. But for temperatures around T ^ 314 MeV in the 
chirally symmetric phase the scalar mixing angle drops 
down to 6*5 ^ —54° (0s ^ 0°). Around these tempera- 
ture the masses of the physical u and /o anticross and the 
ones of the nonstrange-strange (TNS-fs system cross. This 
is displayed in Fig. [T] Hence, for larger temperatures, 
T > 314 MeV the /o is now a purely nonstrange quark 
system and the a a purely strange state. For very large 
temperatures, around 900 MeV, the scalar mixing angle 
turns back to the ideal Bs ~ +35° again and a crossing 
and anticrossing of the corresponding masses takes place 
again. Without anomaly the same phenomenon happens 
qualitatively around some larger temperatures of the or- 
der of T - 325 MeV. 

For finite quark chemical potential and vanishing tem- 
perature the mixing angle show qualitatively a similar 
behavior. Around 71 ~ 350 MeV the pseudoscalar an- 
gle Bp increases towards the ideal value while without 
anomaly the angle is already ideal. In the scalar sector 
the angle is nearly ideal in the broken phase and drops 
down to Qs ~ —54° around \i ~ 500 MeV where again 
the masses of the physical a and /□ meson anticross. 

A finite temperature study of the 77-77' complex includ- 
ing the QCD axial anomaly within a Dyson-Schwinger 
approach and a temperature-dependent topological sus- 
ceptibility can also be found in (i?! . IZsj . By means of 
the Wittcn-Veneziano relation the authors studied the 
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interplay between the melting of the topological suscep- 
tibility and the chiral restoration temperature. The au- 
thors find a strong increase of the 77' mass around the chi- 
ral restoration temperature which makes the extension of 
the Witten-Veneziano relation to finite temperature ques- 
tionable. In the present work a constant anomaly param- 
eter has been used corresponding to a constant topologi- 
cal susceptibility. This means that the U{1)a symmetry 
is not restored around the chiral critical temperature but 
at higher temperatures. 



VI. PHASE DIAGRAM AND THE CHIRAL 
CRITICAL SURFACE 



The phase diagram is constructed in the following way: 
for realistic pion and kaon masses (the so-called physi- 
cal point) the light condensate melts always faster with 
T and/or /x than the strange condensate because the 
strange quark mass is heavier than the light quark mass. 
As the chiral phase boundary we use the inflection point 
in the light condensate. 

Later on, we will also vary the meson masses and cal- 
culate the corresponding phase diagrams. As a conse- 
quence, the ordering of the light and strange condensates 
can be inverted since the kaon mass can become lighter 
than the pion mass. In such cases the strange conden- 
sate drops faster than the light condensate and the chiral 
phase transition is triggered by the strange condensate. 
This has to be taken into account systematically, in par- 
ticular, for the calculation of the chiral critical surface. 
The faster melting condensate has been used in order to 
localize the phase boundary. 
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FIG. 8: Phase diagrams with (7(1) a symmetry breaking for 
different values of m„ = 600 MeV (lower line), 800 MeV, 
900 MeV (upper line). 

In Fig. [8] the phase diagrams in the (T, /x)-plane with 
explicit U{1)a symmetry breaking for three different val- 
ues of rria are shown (lower lines correspond to — 600 
MeV, next lines to nia- — 800 MeV and upper lines to 



TO<T = 900 MeV). For each value of Too- the remaining 
parameters of the model are fitted to the vacuum as de- 
scribed in Sec. IIVI and are kept constant. Since the ex- 
plicit U(1)a symmetry breaking leads only to small mod- 
ifications of the phase boundaries (cf. Figs. [1] and [2|), we 
refrain from presenting the phase diagrams in the absence 
of an explicit U{1)a symmetry breaking term. 

At zero chemical potential a crossover is always found 
due to the explicit symmetry breaking terms. The 
crossover temperature depends on the choice of the sigma 
meson mass. For increasing to^- the pseudocritical tem- 
perature also increases (e.g. for to^- = 600, 800, 900 MeV 
we found a Tc ~ 146, 184, 207 MeV). 

Recent lattice simulations at /i = for three quark 
flavors have obtained values for the pseudocritical tem- 
peratures in the range of Tc = 151(3)(3) MeV [13, 
and Tc = 192(7)(4) MeV These results, apphed 

to the current study, would suggest values for the sigma 
mass in between 600 and 800 MeV. 

At zero temperature a first-order phase transition is 
obtained (cf. Fig. For increasing temperatures the 
first-order transition becomes weaker and terminates in 
the critical end point (CEP). How to measure this point 
and what the distinctive signatures of this point are is not 
yet settled. It is interesting to see that the mass of the 
a meson as a function of temperature and/or chemical 
potential through the CEP always drops below the mass 
of the pion not only for the corresponding two flavor but 
also for the three flavor calculation [2^ . This is a general 
feature of the LaM since the potential flattens at this 
point in radial a direction. In a similar NJL calculation 
this is not the case [i^, [5l| . 

For Too = 600 MeV the location of the CEP is at 
(Tc^c) = (91,221) MeV and for m„ = 800 MeV at 
(63, 327) MeV. As a consequence of the model parameters 
dependence, the location of the CEP moves for increas- 
ing Too towards the /i axis. It is interesting to observe 
that already for m„ = 900 MeV the existence of the CEP 
disappears and the phase transition is a smooth crossover 
over the entire phase diagram. 

Almost no influence of the axial anomaly on the phase 
boundaries and thus on the location of the CEP is seen. 
For comparison, with and without the U{1)a symmetry 
breaking and each for TOo = 600 MeV the location of the 
CEP changes from (r^Mc) = (91,221) MeV to (89,228) 
MeV. For TOo = 800 MeV the changes are even smaller 
(Tc^c) = (63,327) MeV to (63,328) MeV. 

In Ref. [131 a gauged linear sigma model with chiral 
U{Nf) X U{Nf) symmetry without quarks within the 
2PI resummation scheme has been considered. For the 
two flavor case and neglecting the influence of the vector 
mesons the opposite behavior of the chiral phase transi- 
tion as a function of TOo is observed: for /i = a crossover 
is seen for a small a mass (too = 441 MeV) and a first- 
order transition for a large a mass (too = 1370 MeV). If 
the vector mesons are included the transition leads to a 
more rapid crossover and brings one closer to the second- 
order critical point. The conclusion is that the critical 
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endpoint moves closer to the temperature axis. Thus, the 
inclusion of vector mesons should improve the agreement 
with lattice QCD results since usually, in chiral models 
the critical endpoint is located at smaller temperatures 
and larger chemical potentials. 

For three quark flavors renormalization-group argu- 
ments predict a first-order transition in the chiral limit 
independent of the U{1)a symmetry breaking [ssj . This 
behavior is displayed in Fig. [9] where the resulting phase 
diagrams including the anomaly for varying pion and 
kaon masses are shown for rric = 800 MeV. We have 
chosen a path in the (m7r, )-plane through the phys- 
ical mass point towards the chiral limit by varying the 
pion mass and accordingly the kaon mass by keeping the 
ratio mTr/rriK fixed at the value given at the physical 
point m*/m|f. 
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FIG. 9: Phase diagrams witfi U{1)a symmetry breaking 
for rricr = 800 MeV and different pion masses: m-n/m^ = 
0.488 (lower line), 0.6,0.8,1.0,1.2,1.36 (upper line), ml = 
138 MeV, m*K = 496 MeV. The ratio niTr/mK = m^/mjf is 
kept fixed. 

The pion and kaon masses are varied by changing the 
explicit symmetry breaking parameters, and hy, while 
keeping all remaining model parameters fixed at the val- 
ues obtained at the physical point. For a pion mass 
1.36 times larger than the physical one (for mc, = 800 
MeV) the CEP lies exactly on the /x-axis, hence for pion 
masses above this value the phase transition is a smooth 
crossover over the entire phase diagram and no CEP ex- 
ists anymore. 

For a decreasing pion mass, the location of the CEP 
moves towards the T-axis and already for half of the phys- 
ical pion mass the CEP hits the T-axis a.t fi — 0. Re- 
markably also for half of the physical pion mass the CEP 
hits the T-axis in a similar calculation within a two-flavor 
LctM within the same approximation scheme but without 
U{1)a symmetry breaking and for Wcr = 600 MeV [2^ . 
For smaller pion masses the chiral phase transition turns 
into a first-order one for all densities and consequently 
no CEP exists any longer. 

The various transition lines in the phase diagram 



shrink towards the origin of the phase diagram for smaller 
pion and kaon masses because the condensates decrease 
more rapidly as a function of the temperature and chem- 
ical potentials. If one rescales the temperature with the 
critical temperature at /i = and the chemical potential 
with the critical chemical potential at T = then all the 
transition lines in the phase diagram lie on top of each 
other for different pion masses. 

In connection with the existence of the CEP in the 
phase diagram it is interesting to analyze its mass sen- 
sitivity. For this purpose the critical surface for the chi- 
ral phase transition is evaluated in Fig. [10] as a func- 
tion of the pion and the kaon masses with (left panel) 
and without (right panel) the U{1)a symmetry break- 
ing. Again the masses are tuned by variation of only 
the explicit symmetry breaking parameters (see Eq. ([19])) 
while keeping all other model parameters fixed similar to 
Ref. [S^l- For small kaon masses but large pion masses 
the explicit symmetry breaking parameter for the strange 
direction hy in Eq. (|19p can become negative. The cor- 
responding kaon mass, where this happens, is given by 
rriK = "mTT ^/JtiJ^Jk and is shown in both panels as a 
solid line. 

The critical surface is defined by the value of the criti- 
cal chemical potential /ic of the CEP for a given mass 
pair (m^,TOif). It is the surface of the second-order 
phase transition points displayed in a three-dimensional 
(/ic, TOtt, mi<-)-space. For values of the chemical potential 
above the chiral phase transition is of first-order while 
for values below the surface the transition lies in the 
crossover region. With or without anomaly the critical 
surface grows out perpendicular from the (m^, m/^)-mass 
plane at = 0. The tangent plane to the critical sur- 
face has a decreasing slope for larger masses as expected 
from Fig. [5] or Fig. [S] Thus, this model study excludes 
the so-called nonstandard scenario, found in a recent lat- 
tice analysis with imaginary chemical potentials, where 
the first-order region shrinks as the chemical potential is 
turned on [l^ . In the nonstandard scenario the bend- 
ing of the critical surface has the opposite sign and the 
physical realistic mass point remains in the crossover re- 
gion for any fi. In the Figure the physical mass point is 
denoted by an arrow. 

Since the critical chemical potential /ic cannot grow 
arbitrarily the surface must have a boundary and hence 
stops to exist for larger (m^r, T?ij(-)-masses which arc not 
shown in the Figure. 

In Ref. [3^ the LcrM with quarks in an one-loop ap- 
proximation, based on optimized perturbation theory, 
was evaluated and an bending of the critical surface away 
from the m/^-axis at rn^ = for a kaon masses greater 
than 400 MeV was observed. The precise value of tuk 
for the onset of the bending depends on the order of the 
used chiral perturbation theory (xPT) for the baryon 
mass extrapolations. Thus, the unphysical bending indi- 
cates that the validity range of the xPT for the baryons, 
used for the model parameter extrapolation away from 
the physical point, was exceeded. As a consequence, no 
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(a) with U{1)a anomaly (b) without U{1)a anomaly 

FIG. 10: The chiral critical surface in the {mTr,mK) plane for rricr = 800 MeV. The arrow points to the critical quark chemical 
potential at realistic pion and kaon masses (physical point). The solid line in the right panel is given by rriK ~ m,r Ay/ir/2/K. 



(tri)critical point on the m^^-axis for fi = where the 
boundary of the first-order transition region terminates, 
can be located. 

In this work no strong bending of the surface away 
from the m^^-axis is seen for kaon masses not larger than 
600 MeV but it also seems that it never approaches the 
mif-axis at least for mg- = 800 MeV. On the other hand, 
including the anomaly we find a critical m^'''* ~ 177 MeV 
where the critical surface intersects the solid line in the 
left panel of Fig.[Tni On this line and for > m'^'^* the 
transition turns into a smooth crossover. On the con- 
trary, without anomaly (right panel in the Figure), no 
critical pion mass is found at yU = at least for values 
below 200 MeV. This means that the phase transition 
on the solid line is still of first-order similar to the find- 
ings of [U where the infiuence of the anomaly on the 
phase transition for vanishing chemical potentials within 
a SU{3) X SU{3) LaM without fermions in Hartree ap- 
proximation has been investigated. 

Furthermore, the effect of the U{1)a anomaly on the 
shape of the surface is rather marginal for a kaon mass 
greater than 400 MeV (cf. both panels in Fig. [TOl). For 
larger kaon masses the strange sector decouples from the 
light sector and the chiral phase transition is basically 
driven by the (light) nonstrange particles. 

On the other hand, for a kaon masses smaller than 400 
MeV we see a considerable influence of the anomaly on 
the shape of the critical surface. Without anomaly the 
region of first-order phase transitions at ^ = is reduced 
which is in contrast to the results obtained with a LctM 
without quarks js^]. In this reference, it is found that 
the first-order transition region at = grows with and 
without anomaly for increasing sigma masses. For sigma 
masses greater than 600 MeV and without anomaly the 
physical point is well located within the first-order region 
while it is always in the crossover region with anomaly. 

Including quarks we obtain an opposite tendency: the 
physical point is always in the crossover region and for 



larger sigma masses the size of the first-order transition 
region at = decreases as can be seen from Fig. 1111 In 
this figure five crosssections, for rria = 500 . . . 900 MeV of 
the chiral critical surface with (left) and without anomaly 
(right panel) are shown as a function of the pion mass. 
As a trajectory in the (m,r) )-plane a path through 
the physical point towards the chiral limit has been cho- 
sen for these figures. This path is given by fixing the 
pion over kaon mass ratio to the physical one, i.e. set- 
ting m^f/m^ = m*j^/m* where the star denotes the cor- 
responding physical masses. 

For larger ma- values the chiral critical surface with 
anomaly moves to smaller pion masses. This effect is 
less pronounced if the anomaly is neglected (right panel). 
In accordance with Fig. [5] the chiral critical surface line 
for TTicr = 900 MeV terminates before the physical point 
(m,r/m* = 1). Due to the positive parameter we 
cannot evaluate the chiral critical surface for arbitrary 
pion and kaon masses for smaller values of as already 
mentioned. Nevertheless, the results for smaller sigma 
masses are shown as dashed light curves in both panels 
and stop at certain pion mass ratios. 

It is instructive to replace in the last Fig. [TT] the criti- 
cal chemical potential with the critical temperature and 
investigate its mass sensitivity. This leads to Fig. [T^ 
where the critical temperatures of the CEP's as a func- 
tion of the pion mass similar to Fig. [Tl] are plotted for 
five different a meson masses with (left panel) and with- 
out anomaly (right panel). In general, the curves start 
at zero temperature and grow to a certain finite value for 
decreasing pion masses. For the parameter sets without 
spontaneous symmetry breaking the curves with anomaly 
for TTicr = 500 and 600 MeV, and without anomaly for 
rua = 500 MeV, show the opposite behavior. In this 
case, i.e. for smaller values of ma-, and large unphysi- 
cal pion masses mT^/m^ > 1, which is realized in lattice 
simulations, the location of the CEP moves towards the 
T-axis. For decreasing pion masses the critical tempera- 
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FIG. 11: Five cross sections of the chiral critical surface (left panel with, right panel without anomaly) for m^y — 500, . . .900 
MeV. As trajectory in the (m,r,niK) plane we have chosen the path through the physical point towards the chiral limit, i.e. we 
have kept the ratio of 77i,r/mx = m*/m^ fixed, m* and denote the physical masses. See text for details. 





250 r 




200 - 


> 


150 - 










H 


100 




50 








\ 

\ \ 

w 



m„= 500MeV 
m„= eOOMeV 

m„= 700MeV 

m„= SOOMeV 

m„=900MeV 



0.5 



1.5 



> 



250 p 

200 - 

150 - 

100 

50 

- 




m„= 500MeV 

m„= 600MeV 

m„= 700MeV 

m„= SOOMeV 

m„=900MeV 



0.5 



1.5 



(a) with U{1)a anomaly 



(b) without U{1)a anomaly 



FIG. 12: Similar to Fig. 1111 but the critical chemical potentials are replaced by the corresponding critical temperatures. 



tures decrease meaning that the CEP moves towards the 
/Lt-axis in the chiral hmit. Consequently, for these param- 
eter sets no first-order transition occurs in the chiral limit 
as already mentioned. Nevertheless, all of curves inter- 
sect roughly at the physical point (m^/m* = 1). Only 
the extrapolation towards the chiral limit is questionable 
as can be seen by the positive slope of these curves for 
decreasing pion masses. 

VII. SUMMARY 

In the present work chiral symmetry restoration in hot 
and dense hadronic matter is analyzed. As an effective 
realization of the spontaneous breaking of chiral symme- 
try in the vacuum, the SU{3) l x SU (3)_r symmetric LctM 
with quarks has been used. Within this model, the grand 
thermodynamic potential was evaluated in the mean-field 
approximation. Six of the seven model parameters are 



fixed to the low-lying pseudoscalar mass spectrum, which 
is experimentally well-known and to the scalar sigma me- 
son mass. Since the experimental situation in the scalar 
sector is not very certain we have varied the value of the 
sigma mass over a wide range and have investigated its 
consequences for the physics. The remaining model pa- 
rameter, the Yukawa coupling, is determined by the non- 
strange constituent quark mass while the condensates are 
governed by the PCAC relation. 

At the physical mass point for various sigma masses a 
smooth finite temperature chiral crossover at zero den- 
sity and a first-order transition for finite chemical poten- 
tial at zero temperature is found. The U{1)a anomaly 
has only little influence on the strange condensate while 
no modification of the light condensate at the physical 
point is observed. The pseudocritical crossover tempera- 
ture depends on the choice of the sigma meson mass and 
coincides with recent lattice simulations for mass values 
in between = 600 . . . 800 McV. 
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In our approximation no negative squared meson 
masses are generated in the medium as is the case for 
the LcrM without quarks. Low-energy theorems hke the 
Goldstone Theorem or the Ward identities are not only 
fulfilled in the vacuum but also in the medium. This 
enables a careful analysis of the chiral symmetry restora- 
tion pattern of the meson nonets with and without axial 
anomaly. An anticrossing of the scalar cr-/o(1370) masses 
at higher temperature is seen which is reflected in the 
corresponding mixing angle investigation. 

For a sigma masses below 900 MeV a CEP is found in 
the phase diagram. In contrast to a similar calculation 
in the NJL model at this point the sigma meson mass 
drops below the pion mass. This behavior is analogous 
to the findings in the two flavor case. 

The chiral critical surface in the (m7r, ?nif)-planc al- 
ways has a positive slope such that the so-called non- 
standard scenario can be excluded. In the chiral limit 
the expected first-order transition at ^ = is foimd. 
A large modification caused by the anomaly is visible. 
With anomaly the first-order region in the fi = plane 
is enlarged for smaller masses. For larger kaon masses 
the shape of the critical surface becomes independent of 
the anomaly. Furthermore, for increasing ma- the surface 
becomes steeper. 

One drawback of the effective model used here is the 
lack of confinement properties which certainly modify the 
thermodynamics in the broken phase. A first step to- 
wards an implementation of giuonic degrees of freedom 
which could mimic certain confinement aspects can be 
achieved by adding the Polyakov loop to the quark dy- 
namics. Work in this direction is in progress [55|. 
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therefore have used different input values for rria- in the 
range of 400— 1000 MeV and have refitted the remaining 
parameters of the model accordingly. It is remarkable 
that larger sigma meson masses are not adjustable since 
THjfj as a function of the quartic coupling Ai saturates 
around 1100 MeV. 

In Tab. |T] all resulting meson masses are listed. The 
upper block contains the fit without anomaly and the 
lower block the fit including the anomaly. Except for 
the pseudoscalar masses the first four/five columns in 
the table respectively, all other (scalar) masses and the 
mixing angles are predictions of the model. 

In Tab. |TT] the values for the six mesonic model pa- 
rameters of the LcrM with and without explicit U{1)a 
symmetry breaking are summarized. The Yukawa cou- 
pling is always kept fixed to g ^ 6.5 corresponding to 
a constituent quark mass of = 300 MeV. For this 
Yukawa coupling the strange constituent quark mass is 
predicted to be iris ~ 433 MeV. The decay constants, 
92.4 MeV and 113 MeV, are also kept con- 

stant for all fits. It is interesting to realize that for small 
values of rria and with U{1)a symmetry breaking the 
mass parameter changes sign and becomes positive 
when fitted to realistic masses. As a consequence spon- 
taneous symmetry breaking is lost in the chiral limit and 
all condensates will vanish in this limit. This happens for 
77icr < 700 MeV. Even without anomaly breaking (c = 0) 
a similar phenomenon can be seen. For this case the 
masses are smaller when this case sets in, i.e. rria < 500 
MeV. This is the motivation for our choice of nia = 800 
MeV. For this parameter set we can investigate the mass 
sensitivity of the chiral phase transition over arbitrarily 
explicit symmetry breaking values including the chiral 
limit. The choice rric = 800 MeV for the parameter fit 
without anomaly is also in agreement with (33 | where 
Too- — 600 MeV is a misprint in this reference. For larger 
TOct values the quartic coupling Ai increases significantly 
which restricts the parameter fits for larger than 1000 
MeV. 



APPENDIX B: MESON MASSES 



APPENDIX A: PARAMETER FITS 

In this appendix several used parameter fits for the 
linear sigma model (LcrM) with three quark fiavors arc 
collected. 

As experimental input wc have chosen the low-lying 
pseudoscalar meson mass spectrum (tojt, mx, "rriri), the 
constituent quark mass mq and the pion and kaon decay 
constant. To be more precise, for the fit with anomaly 
(c 7^ 0) the sum m^+m^, is chosen as input while for c = 
the TO^ has not been used as input. The experimental 
values, taken from Ref. [i^, are listed in the last line 
of Tab. U for comparison. Since the chiral cr-particle is 
a broad resonance its mass is not known precisely. We 



The scalar = O"*" and pseudoscalar = 0~ meson 
masses are defined by the second derivative w.r.t. the cor- 
responding scalar and pseudoscalar fields tps.a = and 
'Pp,a — T^a (a = 0, . . . , 8) of the grand potential Q.{T,iJLf ), 
Eq. ([7]), evaluated at its minimum. The minimum is 
given by vanishing expectation values of all scalar and 
pseudoscalar fields but only two of them, Ux and ct^, are 
nonzero. 



s,p 



(Bl) 



In the vacuum the contribution from the quark poten- 
tial vanishes. Hence, only the mesonic part of the poten- 
tial determines the mass matrix completely. The squared 
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niK 






9p 


mag 






Os 


400 


138 


496 


138.0 


634.8 


35.3 


850.4 


1124.3 


1257.3 


16.7 


500 


138 


496 


138.0 


634.8 


35.3 


850.4 


1124.3 


1267.4 


18.7 


600 


138 


496 


138.0 


634.8 


35.3 


850.4 


1124.3 


1282.3 


21.5 


700 


138 


496 


138.0 


634.8 


35.3 


850.4 


1124.3 


1304.9 


25.5 


800 


138 


496 


138.0 


634.8 


35.3 


850.4 


1124.3 


1341.4 


31.3 


900 


138 


496 


138.0 


634.8 


35.3 


850.4 


1124.3 


1408.0 


40.0 


1000 


138 


496 


138.0 


634.8 


35.3 


850.4 


1124.3 


1563.4 


53.2 


400 


138 


496 


963.0 


539.0 


-5.0 


1028.7 


1124.3 


1198.4 


14.9 


500 


138 


496 


963.0 


539.0 


-5.0 


1028.7 


1124.3 


1207.5 


16.9 


600 


138 


496 


963.0 


539.0 


-5.0 


1028.7 


1124.3 


1221.1 


19.9 


700 


138 


496 


963.0 


539.0 


-5.0 


1028.7 


1124.3 


1242.3 


24.2 


800 


138 


496 


963.0 


539.0 


—5.0 


1028.7 


1124.3 


1278.0 


30.7 


900 


138 


496 


963.0 


539.0 


-5.0 


1028.7 


1124.3 


1348.0 


40.9 


1000 


138 


496 


963.0 


539.0 


-5.0 


1028.7 


1124.3 


1545.6 


57.1 


400-1200 


138.0 


496 


957.78 


547.5 




984.7 


1414 


1200-1500 





TABLE I: Meson masses and mixing angles in the vacuum for different sets of parameters. The first six columns are input 
while the remaining columns are predictions. Upper block: without (c = 0) U{1)a anomaly, lower block: with U{1)a anomaly. 
Last line: experimental values from the PDG [42| . 



ma[MeV] 


c[MeV] 


Ai 


m2[MeV^] 


Aa 






400 





-24.55 


-f(309.41)2 


82.47 


(120.73)3 


(336.41)3 


500 





-21.24 


+ (194.82)2 


82.47 


(120.73)3 


(336.41)3 


600 





-17.01 


-(189.85)2 


82.47 


(120.73)3 


(336.41)3 


700 





-11.61 


-(360.91)2 


82.47 


(120.73)3 


(336.41)3 


800 





-4.55 


-(503.55)2 


82.47 


(120.73)3 


(336.41)3 


900 





5.56 


-(655.82)2 


82.47 


(120.73)3 


(336.41)3 


1000 





24.22 


-(869.50)2 


82.47 


(120.73)3 


(336.41)3 


400 


4807.84 


-5.90 


-f(494.75)2 


46.48 


(120.73)3 


(336.41)3 


500 


4807.84 


-2.70 


+(434.56)2 


46.48 


(120.73)3 


(336.41)3 


600 


4807.84 


1.40 


+(342.52)2 


46.48 


(120.73)3 


(336.41)3 


700 


4807.84 


6.62 


+(161.98)2 


46.48 


(120.73)3 


(336.41)3 


800 


4807.84 


13.49 


-(306.26)2 


46.48 


(120.73)3 


(336.41)3 


900 


4807.84 


23.65 


-(520.80)2 


46.48 


(120.73)3 


(336.41)3 


1000 


4807.84 


45.43 


-(807.16)2 


46.48 


(120.73)3 


(336.41)3 



TABLE II: Different parameter sets for various with (c 7^ 0) and without (c = 0) (7(1)a anomaly. 



mass matrix is diagonal and due to isospin SU(2) sym- 
metry several matrix entries are degenerate. We begin 
with the scalar, = 0"*", sector corresponding to z = ,s. 
The squared mass of the ao meson is given by the (II) 
element which is degenerate with the (22) and (33) ele- 
ments. Similar, the squared k meson mass is given by the 



(44) element which is also degenerated with the (55) , (66) 
and (77) elements. The a and /o(1370) meson masses are 
obtained by diagonalizing the (00)-(88) sector of the mass 
matrix introducing in this way a mixing angle 9s- Ex- 
plicitly, the squared masses for scalar sector, formulated 
in the nonstrange-strange basis, are 
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with 



2 2 

= m 



2 ' 

Ai(a2 + dl) + + ^dxOy + 2aj^ 



c _ 



*s,00 



rn^ 00 cos^ 6*5 + m^ ^g sin^ 6*5 + 2m^_(,8 sin cos 6*5 , 
00 sin^ + cos^ Bs — '^^1,08 sin ^5 cos 9s 

+ y (7^^ + ^V2axay + 5ct^) + A2(ct^ + a^) - ^(%/2ct^ + ay) , 



3V2 



(B2) 
(B3) 
(B4) 
(B5) 



(B6) 



The situation for the pseudoscalar sector (i ~ p) is 
completely analogous with the following labeling: the 
squared pion mass is identified with the (11) element and 
the squared kaon mass with the (44) element of the pseu- 
doscalar mass matrix. Similar to the scalar case, the rj 



and 77' mass are obtained by diagonalizing the (00)- (88) 
sector and accordingly a pseudoscalar mixing angle 0p is 
introduced. Explicitly, the squared masses for the pseu- 
doscalar sector are 



ml 


= + 






U •^2-2 


V2c 
— o-y 


(B7) 


rriK 


= + 




-D- 




- \/2ax&y + 2afj - 


(B8) 






cos^ 9p 


+ ml 


gg sin^ 9p 


+ 2TOp 08 sin^p cos 9p 


(B9) 


m\ 




sin^ 9p 


^ ml^gg cos^ 9p 


— 2inp ogsin^pcosflp 


(BIO) 




= + 




-D- 




' al) + ^{2ax + V2ay) 




'^p,88 


= + 


Ai(a2 + 






^ Ad-D ~ ^{Aa, ~ V2<jy) 





<08 -^i^l- 2^1) ~ ^(V2a. - 2ay) 



Both mixing angles are given by 
2mf Qg 

tan20i — — 2 ' — 2 — ' * ~ ^^P ■ 

"^1,00 ~ "^1,88 

In the medium the meson masses are further modified 
by the quark contribution ([5]). In order to evaluate the 



second derivate (|Bip for the quark contribution the com- 
plete dependence of all scalar and pseudoscalar meson 
fields, cf. Eq. ([3]), in the quark masses has to be taken 
into account. The resulting quark mass matrix can be 
diagonalizcd. Finally, we obtain the expression 



min f=q,s' 



(27r)3 2E,j 



2</ . 



(B12) 
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TTl 
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714 


7r4 
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TTS 


TTS 
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3 
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3 
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TABLE III: Squared quark mass second derivatives with re- 
spect to the meson fields evaluated at the minimum. Left 
column block, the sum over two light quark flavors, denoted 
by index / and right column block only the strange quark 
flavor, index s. 



where 



we have introduced the short 
for the quark 



= drnj/dipi^a 



hand notation 
mass derivative 



w.r.t. the meson fields f iM, the quark function 



(B13) 



and correspondingly the antiquark function bqj{T, fif) = 
bqj{T, —^f). The index i distinguishes between a scalar 
and pseudoscalar field which we omit in the following. In 
Tab, mil all second quark mass derivatives w.r.t. the me- 
son fields replaced by the non- vanishing vacuum expecta- 
tion values in the nonstrange-strange basis are collected. 
Despite the SU{2) isospin symmetry the derivatives are 
different for the up- and down-quark sector. In the table 
mil left column block, the sum over the two light quark 
flavors is shown which leads to large cancellations. 



These states are associated to each other by a rotation 
with an angle a = — arctanv^ ^ —54.74° 



I '7ns; 

\vs) 



V3 ( 



1 V2 
/2 1 



\V8) 

\Vo) 



(C3) 



Diagonalization of the mass matrix in the (Tys,^) ba- 
sis is achieved by the introduction of the pseudoscalar 
mixing angle 9p. This yields the relations 



w 



cos tip — sm f p 
sin 9 p cos 9 p 



\V8) 

\vo) 



(C4) 



For the (?7nSi '7s) basis a similar relation with the mixing 
angle (j)p holds 



\v) 
W) 



cos I 
sin f 



- sm t 
cos t 



l'?Ns; 
\vs) 



(C5) 



For vanishing mixing angle (j)p corresponding to 9p = 
— arctan %/2 -54, 7° the 77 tends to a pure nonstrange 
r^NS Eind rj' to a pure strange rjs. In contrast, for a mixing 
angle (f>p = 90° {9p ~ -1-35, 3°) the ordering is reversed 
and rj r/s and 77' r/NS- The ordering transition 

occurs at (t>p = 45° {9p 9.74°). 

The diagonalization of the mass matrix in the (?7ns7 Vs) 
basis leads to the masses 



- sin^ ( 

TO^g COS^ ( 



)p — m: 

}p + m. 



2 



and to the mixing angle (j)p given by 
tan 2(f>p 



sin2(2,/)p),(C6) 
sm^{2(t)p). (C7) 



(C8) 



ns 



r/NS 



Equivalently, these expressions can be rewritten in a form 
which do not contain the mixing angle explicitly 



APPENDIX C: ISOSCALAR MIXING 

In this appendix a collection of relations describing 
the mixing of isoscalar states in the pseudoscalar and 
scalar multiplet is presented. The isoscalar (/ = 0) pseu- 
doscalar states in the octet-singlet (778,770) basis are de- 
fined by 



= — r= \uu + dd - 



2ss) , \T]n) = |uM - 



dd + ss) . 
(CI) 

For a realistic flavor breaking in the vacuum the physical 
77 meson is close to the 773 and 77' to ryg. 

The eigenstates in the flavor nonstrange-strange 
iv^SiVs) basis are given by 



\VNs) 



1 

71 



uu + dd) , \ris) = \ss) 



(C2) 



1. 



2 

'/NS 



I, 

o('« 



2 



')S 



), 
), 



with 



(C9) 
(CIO) 

Note, that 

these expressions are numerically more stable compared 
to (|C6I - IC7|) because possible ambiguities in the tangent 
()C8|) do not appear here. 

The matrix elements in the (TyNSi^ys) system are ob- 
tained by a base change from the ones in the (r]g,riQ) 
basis with the result 



JJNS 



VS 



) +4?7l2 



ms,vs 



m 



2 



VS 



(2'"p,oo 

('"p.OO -+ 



\{V2{ml, 



>,88 



2V2m^ 



p,08J 1 



2m2 «« - 2\/2m^_o8) , (Cll) 
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As a consequence, the mixing angles (j>p and 9p are 
related by 

(/)p^ep + arctan V2 ^ Op + 54.74° . (C12) 

Furthermore, supposing < m^i one finds with 
(IC6l),(lC7| for < 45° {Op < -9.74°) m^^^ < m^^ 
while for (j)p > 45° the ordering of the masses in the 
nonstrangc-strange system are reversed. 

Scalar mesons differ from the pseudoscalar ones only 
in the orbital excitation. Hence, all quoted relations can 
be immediately converted to the scalar (cr, /o) complex 



with the corresponding replacements, e.g. 

n/o)\ /cos0, -sin^,\ naNs)A 
\\a) J ysm0s cos 0s / V I'^s) / 

For an ideal scalar mixing angle (j>s = 90° the a meson 
is a pure nonstrange state cns and /o —crs- Further- 
more, a matches with rj' and /o with ?/. 

Since the mass of the /o meson is larger than mo- we 
obtain for an ideal mixing (jjs ~ 90° the ordering TOq-s > 
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